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Abstract 


We present a new model for calculating locomotion costs in breath-hold divers. Starting from basic mechanics principles, we 
calculate the work that the diver must provide through propulsion to counterbalance the action of drag, the buoyant force and 
weight during immersion. Compared to those in previous studies, the model presented here accurately analyses breath-hold divers 
which alternate active swimming with prolonged glides during the dive (as is the case in mammals). The energy cost of the dive is 
strongly dependent on these prolonged gliding phases. Here we investigate the length and impacts on energy cost of these glides 
with respect to the diver characteristics, and compare them with those observed in different breath-hold diving species. Taking into 
account the basal metabolic rate and chemical energy to propulsion transformation efficiency, we calculate optimal swim velocity 
and the corresponding total energy cost (including metabolic rate) and compare them with observations. Energy cost is minimised 
when the diver passes through neutral buoyancy conditions during the dive. This generally implies the presence of prolonged 
gliding phases in both ascent and descent, wh ere the buoy a ncy (varying with depth ) is best used against the drag, reducing energy 
cost. This is in agreement with past results ( Miller et al. . 20121 : Sato et all 2013h where, when the buoyant force is considered 
constant during the dive, the energy cost was minimised for neutral buoyancy. In particular, our model confirms the good physical 
adaption of dolphins for diving, compared to other breath-hold diving species which are mostly positively buoyant (penguins for 
example). The presence o f prolonged glides implies a non-t r ivial d ependency of optimal speed on maximal depth of the dive. This 
extends previous findings ( Sato et all 20101 : Watanabe et al. , 2011 ) which found no dependency of optimal speed on dive depth for 
particular conditions. The energy cost of the dive can be further diminished by reducing the volume of gas-filled body parts in 
divers close to neutral buoyancy. This provides a possible additional explanation for the observed exhalation of air before diving in 
phocid seals to minimise dive energy cost. Until now the only explanation for this phenomenon has been a reduction in the risk of 
decompression sickness. 
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1. Introduction 


During their dives, breath-hold diving animals minimise en¬ 
ergetic cost to gain time foraging as oxygen stored in their 
body is limited. Besides plastic physiological adaptations to 
diving, like brady cardia, reduction and redistribution of the 
blood flow, etc., ( Butler and Jonesl 1997: Koovmanl 1985 : 


Koovman and Ponganisl 1998t iButlei . 12004 ). dive energy cost 
can be lowered by reducing dive duration and/or the mechan¬ 
ical work necessary for propulsion. Energy cost related to 
the basal metabolic rate is proportional to dive duration and 
inversely proportional to swimming velocity. On the other 
side, energy spent for propulsion depends on the drag force 
during the dive, which increases with the square of veloc¬ 
ity. Besides swimming optimisation and hydrodynamics, thrust 
work is efficiently reduced by slowing down swim speed. 
Optimal dive velocity is a compromise, taking into account 
these two elements and is specific to the diver’s body char¬ 
acteristics and maximal dive depth. An additional energy 
cost reduction is obtained, particularly in mammals, by al- 
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ternating active swimming with prolonged glides using buoy¬ 
ant force and weight to their advantage but also by varying 
their s troke frequency and/or by using s t roke-and-glide swim- 


ming 

(Crocker et all 

1997; 

Webb et all 1998; Skroyanetal. 

1999; 

Williams et al.. 

200C, 

Dayis et all 200 It Nowacek et al. 

2001, 

Biuw et all 120031: ISato et all 20031: Miller et a 

1 2004 

Watanabe et all 12006; Aoki et al.. 2011; Maresh et all 

2014. 


In recent years many models have been developed to 
study energy ec onomy in bre a th-hol d divers, mainl y birds 


and mammals ( 

Wilson et all 19921 

Skroyan et al.. 

1999 

Hansen and Ricklefs. 

2004; Sato et alE 

lOlOl Aoki et al.. 

2011 

Watanabe et all 2011 

; Miller et all 2012b. 1. In these stud- 


ies, energy cost is estimated from first principles from the 
work of the mechanical forces acting on the diyer and the 
diyer’s metabolism. A complete force and energy cost anal¬ 
ysis has been made for penguins, which are usually posi- 
tiyely buoyant and glide only during the ascending phase. 


(Hansen and Ricklefs 

. 2004; Sato et al. 

. 2 OIOI Watanabe et al. 

201 ll iMillere 

tal., 2012h. In particular Sato et al. (l2010b: 

Watanabe et al. 

(I 2 OII) found the analytical expression of op- 


timal yelocity for minimising energetic cost of the diye for this 
particular case. 
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In the case of mammals, the problem is more complex due 
to their alternation of prolonged glides and active swimming 
phases - crucial for energy balance. For this case, an exhaus¬ 
tive analy sis is still missing in the literature. An investigation 
similar to Sato.eL^ (120101) for mammals was developed by 
Miller et alJ ( 2012*1) 1 though the de pendency of buoyan cy on 
depth was not taken into account. ISkrovan et al.l ( 1999h pro¬ 
duced another analysis where only the average value of the 
buoyant and thrust forc es was considered . A nu merical study 
was also developed by Davis and WeihsI (l2007h for elephant 
seals, but in this case it was for a related investigation where 
the benefit of transit dives with gliding descent is compared to 
purely horizontal displacement. 

Here we present a new model for breath-hold diving for a 
general case which can include the presence of prolonged glid¬ 
ing phases in both ascent and descent. We begin with analysis 
of the mechanical forces involved and their exact dependency 
on depth. We calculate length of the gliding phases and take 
it into account to evaluate dive locomotion cost. Moreover, by 
including some basic consideration of the animal physiology, 
we also predict the energy cost of a typical dive. We study 
the dependency of locomotion and energy cost on the differ¬ 
ent parameters (dive velocity, diver buoyancy, mass, etc.) and 
their optimisation. We then compare this with observations of 
real dives in different breath-hold diving animals, in particular 
the bottlenose dolphin, for which several vertical-dive measure¬ 
ment from trained individuals are available 

In the following section of this paper we present our model 
to calculate mechanical work and total energy cost required for 
a typical dive as a function of the diver and the dive charac¬ 
teristics, including animal metabolism. In the third section we 
discuss our results and how they compare with observations of 
real dives. Section four is our conclusion. 


2. Model description 


2.1. Basic assumptions 

We calculate dive energy cost from first principles from 
the work of the mechanical forces acting on t he diver 


Skrovan et alJ ( 1999b; 

Hansen and Ricklefs 


2004); 

Sato et al. 

(201C 

); Aoki et al.l (l201 1); 

Watanabe et al. ( 

2011); 

Miller et al. 


(120121) . During a dive in a fluid (sea water in our case), the an¬ 


imal is subject to four forces; drag, weight, the buoyant force 
and thrust (see scheme in Figure [TJ. Drag is caused by fric¬ 
tion and the pressure gradient differential. It is strongly de¬ 
pendent on the body velocity v = |v| and always opposes it 
(Fd = -Fdiv)v with V = v/v). 

Weight is proportional to body mass m and the acceleration 
of free fall g {Fy, - mg). The buoyant force is its counter¬ 
part and is proportional to body volume V and the fluid density 
p {Fi, - -pVg). The buoyant force changes with depth due 
to the presence of gas-filled parts of the diver’s body (lungs, 
oral cavities, etc.) whose volume Vg varies with the pressure P 
{PVg = const.). To take into account the effect of Vg variations. 



Figure 1: Scheme of the forces acting on a breath-hold diver (a dolphin in the 
sea in this example) moving with a velocity v in a liquid. Forces are: thrust Ft, 
drag Fd, buoyant force F^ and weight Fw Vg indicates the volume of the gas- 
filled parts of the diver body (mainly the lung, oral cavities, etc.) compressible 
with pressure changes. Vt indicates the volume of the liquid and solid parts 
of the diver body (body tissues), which do not change significantly in volume 
when the pressure increases. 


it is better to rewrite Ft, as a function of the depth d as: 


Fb ^p(Vt + Vg)g = mg 


R,- 


R. 


1 H- ^ 

Po 


( 1 ) 


where V, represents the volume of the tissue of the diver’s body 
comprising liquids and solids which is considered constant at 
any depth, due to its small variation with pressure, as well as 
the water density p. We also introduce the ratios R, = pVtIm 
and Rg - pVgIm, where Vg represents the gas-filled body parts 
volume at the surface. /?, can also be written with respect to 
fluid density p and average body tissue density p,: Rt - pjp,. 
The use of the ratios R, and Rg instead of the volumes V, and 
Vg has the advantage of allowing us to investigate the depen¬ 
dency of locomotion cost and other pertinent quantities on the 
different parameters independently (or almost) from the specific 
mass of the diver. 

Locomotion cost required to descend to a maximal depth D 
and return to the surface can be quantified by the calculation 
of the work W of the thrust force produced by the diver along 
to the dive path. In order to reduce W, breath-holding animals 
(and also human freediver athletes) use the buoyant force and 
weight in their favour by gliding for as long as possible. In 
particular, some mammals such as dolphins, seals, etc. glide for 
large parts of both descending and ascending dive phases. For 
these mammals, a typical dive is composed by a first step where 
the diver makes an effort to descend from the surface. Then 
it actively swims down to a depth do (descent critical depth), 
where the buoyancy is small enough to allowing it to glide to 
the maximum depth D. Following this, the diver makes another 
effort to ascend from the maximum depth to a second depth dA 
(ascent critical depth) from which it glides to the surface under 
the action of the buoyant force. The mechanical work for the 
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Symbol 

Meaning 

Units 

where d( is the infinitesimal of the dive path. 

m 

mass of the diver 

kg 


D 

maximum depth of the considered 

m 

2.2. Additional assumptions and approximations 

T 

total diving time 

s 

Several simplifications and assumptions are necessary to 

V 

swim velocity 

m s“' 

study this formula and compare its predictions with observa- 

Fd 

drag force 

N 

tions. These assumptions are discussed in the following sec- 


weight 

N 

tions. 

Fb 

buoyant force 

N 


Ft 

thrust 

N 

2.2.7. Trajectory and swim velocity during the dive 

C 

drag constant 

Kg m“' 

To calculate dive locomotion cost, we consider purely ver- 

A 

ratio between active to 
passive swimming drag 


tical dives, without horizontal displacement and without con¬ 
sidering possible foraging at the bottom or intermediate depths. 

g 

free fall acceleration 

m s“^ 

We assume in addition that the diver travels at a constant ve- 

Po 

pressure at the surface 

Pa 

locity for the descent and the ascent, which is partially justi- 

P 

liquid density (here the sea density 
assumed to be 1027 Kg m“^) 

Kg m“^ 

fied by the observation of velocity profiles of many breath-hold 
mammals’ diyes (Skroyan et al.. 1999; Williams et al.l 2000t 

Pt 

diver body tissue average density 

Kg m“^ 

Biuw et al.l 2003). We must note that diyes in natural con- 


volume of the body parts 
filled with gas (compressible) 

m3 

ditions are rarely exclusiyely yertical and they may also haye 
a non-constant swim yelocity profile. But for a giyen maxi- 

yo 
^ g 

volume of the gas-filled 
body parts at the surface 

m3 

mal depth, a yertical trajectory minimises the locomotion cost. 
Analysis of yertical diyes then proyides maximal or minimal 

Fs 

= pVgfm, ratio between the 
gas-filled body volume and the mass 


yalues of the inyolyed quantities (total locomotion cost, gliding 
phase extension, etc.) for more general diye conditions. More- 

V, 

volume of the body tissue 
(assumed incompressible) 

m3 

oyer, understanding this simple case proyides the basis for a 
possible generalisation to diyes with an arbitrary diye pitch an- 

Rt 

= pVtfm - pipt, ratio between the 
tissue body volume and the mass 


gle and yelocity profile. 

The assumption of a constant yelocity implies no accelera- 

do 

depth where Ft + Fd - F„ 

(descent critical depth) 

m 

tion, i.e. the resulting total force is equal to zero. Giyen that, 
the thrust Ft must balance the other forces with 

dA 

depth where F„ + Fd - Fb 
(ascent critical depth) 

m 

Ft(cT} = -(Fb{d) + F„ + Fd), (3) 

Ik 

work of the thrust for the dive 
(locomotion cost of the dive) 

j 

where dependency on depth comes only from buoyancy. 

Sp 

efficiency to transform muscular 
movement into propulsion 


Iln the gliding phases of the diye, yelocity tends to increase 
due the effect of buoyant force yariation with depth. It could be 


efficiency to transform chemical 
energy into muscular movement 


kept constant with an intentional increasing of drag without any 
additional significant effort. A deyiation from the assumption 

S 

= SpSm, efficiency to transform 
chemical energy into thrust 


of constant yelocity during the glides does not, howeyer, affect 
the consideration of locomotion cost. On the contrary, it can 

E 

= Em + E,, total chemical energy 
required for the dive 

j 

influence the total metabolic energy cost eyaluation for the diye, 
which could be oyerestimated. 

E, 

chemical energy required 
for the thrust 

j 

2.2.2. Swim technique and accelerations 

E,n 

chemical energy required for 
the basal metabolism 

J 

Besides actiye swimming with a continuous series of strokes 
or prolonged glides, a third swim mode, called stroke-and-glide 

B 

basal metabolic rate 

w 

(a short acceleration phase produced by a single or a small se- 

COT 

= EjirnD), cost of transport 

3 

1 

OQ 

1 

ries of strokes followed by a small glide), is commonly ob- 


Table 1: List of symbols with corresponding units. Bold symbols indicates vec¬ 
torial quantities. The corresponding symbols not in bold indicate their absolute 
values. All quantities are in the units of the international system of units (SI). 
Quantities without units are dimensionless. 


round-trip can be formally written as 

WiD)^ f°Ft(d) de+ Ft(d)-de, (2) 
Jo Jd 


seryed in breath-hold diyers (Webb et al., 1998; 

Skroyan et al. 

1999 

; Nowacek et al. 200 

; Sato et all 2003 

; iMiller et al. 

2004 

; Watanuki et al.. 20051 

Watanabe et all 20061 Aoki et al. 

2011 

; Satoetal.. 2011; Shiomietal.. 2012). 

Here we con- 


sider only the resulting averaged thrust and the corresponding 
average velocity neglecting possible short local accelerations 
present in the stroke-and-glide swim mode. In the same way, 
we do not consider the work required at the beginning of the 
dive to accelerate to the cruise velocity and, at the maximum 
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depth, to change direction, which is about 3/2 mv^, generally 
much smaller than the total work. 


2.2.3. Passive and active drag 

For typical swim velocities and body masses of breath-hold 


diver s, turbulent regime must be considered (IGazzola et al. 
2 OI 4 I) . Drag is then proportional to the square of the veloc¬ 
ity, Fd{v) - Cv^. C is a constant that includes dependency 
on diver body characteristics and the Reynolds number at this 
velocity regime. Possible additional deviations of drag depen¬ 
dency on the square of velocity are not considered. For a given 
cruise velocity, the drag force can be different if the diver is ac¬ 
tively swimming or is gliding. Generally this difference is rep¬ 
resented by the parameter A equal to the ratio of active to pas- 


(Liehthill, 1971; 

Webl 

i 1975: Eish. 

1998: Skrovanetal.. 19991 

Anderson et all 2001 

lEish et all 2014). However, studies that 


include the metabolic rate for locomotion and thermoregula¬ 
tion at the same time indicate that the drag during the active 
swimm ing phases should be sm aller than during the gliding 
phases ( Hind and Gurney . 1997h with 0.2 < d < 1 (depending 


on the diver species). Experiments with a robotic fish mode l 
measure also produce a value of d < 1 ( Barrett et al.L 1999h . 


which is con firmed by recent numerical invest igations on hy¬ 
drodynamic ( Boraziani and Sotironoulosi 2008 ). In our model. 


we c onsider the conserv ative (and simple) assumption d = 1, as 
in to lMiller et al.l j2012h . A deviation from this assumption will 
be discussed in the following sections. 


Unlike [Hansen and Ricklefsl (120041) . which investigated 


Adelie penguin dives, we do not take into account a possible de¬ 
pendency on depth d of the drag constant C (passive to active). 
This dependency is particularly important for divers with a large 
portion of their total volume filled with air or other gases, as in 
the case of penguins. Mammals, the particular focus of this ar¬ 
ticle, have a much smaller relative gas-filled volume, and so are 
much less subject to variations of C with respect to d. 


2.3. Locomotion cost calculation 

From Eqs. dDEl), and the considerations above, the propul¬ 
sion work can be written as 


W(D) ^ [Fd - mg{l - R,)]dD+ 

niRgPo I 
+ ^^log 1 


Pgdo 


+ [Fd + mg{l - R,)](D - dA)+ 
mRnPo 


■log 


ro 


1 -H 

Pb J 


(4) 


where m is the diver mass and Rg and R, are the gas-filled diver 
body volume and tissue volume over mass ratios, respectively. 
Pq is the pressure at the surface, p is the water density, g is 
the free fall acceleration constant and do and dA are the critical 
depths where the diver starts to glide in the descent and the 
ascent, respectively. 


Depending on the maximum depth of the dive D and the diver 
characteristics, the critical depths exist if 

Rg Fd Fd 

Rf H- 77 H- < 1 < Rf + Rp H-for do, (5) 

' 1 + mg ^ mg 

Rg Fd Fd 

R, + - —<1<R,+Rg-— fordA. (6) 

1 + mg ^ mg 

If these conditions are satisfied, the descent and ascent critical 
depths do and dA are 


dD.A - 


pg i-R,^^g 


(7) 


By observing whether or not the critical depths exist, we can 
distinguish different cases of dive: 

Case A Both do and dA exist (the diver partially glides in the 
descent and ascent). 

Case B Only dp exists (the diver never glides in the ascent). 

Case C Only dA exists (the diver never glides in the descent). 

Case D No critical depth exists and F^^, > Ft + Fd for any 
considered depths (the diver is always negatively buoyant, 
and glides for the entire descent). 

Case E No critical depths exists and Ft > F„ + Fd for any 
considered depth (the diver is always positively buoyant, 
and glides for the entire ascent). 

Case F No critical depths exist and Ff > m?c>dFt-F„, F„-Ft) 
for any considered depth (no gliding phase is present, the 
drag is dominant). 

In case A, where both critical depths exist, from the combi¬ 
nation of Eq. (in with Eq. (I?]) we can obtain the explicit formula 
of W. Similar expressions of W for cases B and C can be ob¬ 
tained by substituting one of the critical depths with one of the 
extremes of the [0, D\ interval: dA - 0 for case B and do - D 
for case C. 

Eor cases D and E, the diver glides for the entire descent or 
ascent, respectively, and it never glides in the opposite phase 
{do and dA do not exist) leading to the simple formulae 


W^-^(D) = FdD ± mg(l - R,)D + 


mPoRg 


log 1 


f )■ <*' 


Case E is a peculiar instance where the drag is so strong (or 
the diver is so close to neutral buoyancy and Rg <& 1) that no 
gliding phase is present at all. The work expression for this 
case is even simpler than cases D and E and it can be formally 
obtained from W^’ + leading to 

W^(D) = 2FdD. (9) 


We note that in this case the work for a vertical dive is the same 
that as for an equivalent horizontal displacement, where only 
the drag is relevant. Case E is also complementary to case A. 
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The selection of the case (A or F) that must be considered de¬ 
pends on the values of the various parameters It can be demon¬ 
strated that if 


pgD mgRg 
2Po I+ ef 


( 10 ) 


the existence of both do and dA simultaneously in the interval 
[0, D\ is not possible for any value of Rt and so case A is not 
allowed. 


2.4. Total metabolic energy cost and its optimisation 

To consider the metabolic cost of the dive, rather than lo¬ 
comotion cost, the relevant quantity is total energy E which 
includes basal metabolic rate and the efficiency to transform 
chemical energy to thrust. E can be broken down into two parts 
E - E„, + E,. Em - B T is proportional to basal metabolic 
rate B and the elapsed time during the dive T. E, - W/e is 
proportional to dive mechanical work over the efficiency e to 
transform the chemical energy into forward thrust, s = EpEm 
depends on the propulsion efficiency Sp to transform muscular 
movement into thrust and on metabolic efficiency e,„ to trans¬ 
form chemical energy into muscular work. 

Following these considerations we can write the expression 
for total energy E cost for a dive as: 


2D 


E(D, v)^B — -H 

V 


W{D, v) 

E 


( 11 ) 


W{J) 



W(J) 



remembering that we assume a constant swim velocity v during 
the dive whose resulting duration is T = 2D/v and where we 
consider that efficiency e does not depend on the swim speed. 
In addition, we make explicit the dependency of thrust work W 
on the velocity. 


Figure 2: Thrust work for a bottlenose dolphin of mass m = 177 Kg and drag 
constant C = 4.15 kg diving to a maximum depth of 100 m with a speed 
of 1.8 m s“^ Top: Variation of W with respect to Rf keeping Rg at the constant 
value of 0.0493 (= Bottom: Variation of W with respect to Rg keeping 

Rt at the constant value of 0.9809 (= Different colours con'espond to 

different diving cases. 


3. Results and Discussion 


3.1. Locomotion work 

Some general considerations can be made on the formulae 
of W for the different cases. For any expression of VT a linear 
contribution of Fj is present, as for a horizontal displacement. 
For the simple cases D and E a linear dependency on m, Rg and 
Rt also exists. For cases A, B and C there is a more complex 
dependency on these parameters because of the presence of dojt 
(Eqs. O). Variation of Rg,Rt and EjUmg) values imply the 
existence, or not, of the two critical depths and so determines 
the selection of the case which must be considered for the work 
calculation. 

Compared to the other parameters, only Fj can be easily var¬ 
ied modulating the swim speed. Both Rg and Rf are specific to 
the diver’s body characteristics and have a strong influence on 
the locomotion cost W. The ratio Rg - pVgIm can be modu¬ 
lated by inhaling more o r less air before diving and is gener- 
ally much less than one (IKoovman and Ponganisl Il998h . The 
ratio Rt, equal to the ratio between the sea water and average 
body tissue densities pipt, is generally close to one. Rt can be 
tuned by varying the dive r’s bodily lipid content, as seasonall y 


tuned by varying me dive r s bodily lipid content, as seasonall y 
happens in elephant seals ( Biuw et al.L 2003 : Aoki et al. . 201 ih . 


Erom Eqs. (|7]l we see that small variations of Rt can drastically 
change the value of the critical depths and dive work. 


To compare the model predictions with observations, we con¬ 
sider the case of bottlenose dolphin, for which observations 
from trained individuals p erforming purely vert i cal dives with¬ 


out fo raging are available ( Williams et al. . 1999t Skrovan et al 


1999h . i.e. the ideal conditions in which to test our model, and 


where all cases A-E are seen. To investigate the dependency 
of the dive work on these important parameters Rt and Rg, we 
consider W for a particular bottlenose dolphin individual with 
a mass of 177 Kg, lung volume of 8.5 1 and where we arti¬ 
ficially change the values of Rt and Rg to the values deduced 
from the observation s = 0.9809 and - 0.0493 

( Skrovan et al.[ 1999h . We also ha ve the value of the dra g con¬ 
stant value C = 4.15 kg m“* from ISkrovan et al.l (Il999h for the 
calculation of Ed. 


In Eigure (top) we consider the locomotion cost for a ver¬ 
tical dive to a maximal depth of 100 m, with a swim velocity 
of 1.8 m s ', with Rg - and where we vary Rt. Different 

values of Rt correspond to different cases of dive. Eor this par¬ 
ticular parameter set, case A is possible (both critical depths can 
be present for a specific dive). Similarly, we study the depen- 
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dency of W on Rg varying its value keeping R, - constant 
(Figure|2]bottom). For variations of Rg we do not take into ac¬ 
count its minimal value due to the residual lung volume that, 
due to alveolar collapse, can be very small at high pressure s, 


(IKoovman et allll972l : lFalke et allll985l : lMoore et al.Ll2Qllh . 

In both cases, a minimum of W is present at - 0.9871 
and Rg'^^' - 0.0578 when Rg and R,, respectively, are kept hxed, 
not far from the real values of the considered dolphin. Both 
minima, are included in case A, when both critical depths ex¬ 
ist. More generally, it can be demonstrated that, when case 
A is possible, the minimum of W is unique and is always in¬ 
cluded in case A. The corresponding critical depths are around 
the middle-depth deqOA - Djl + h.deq, where h.deq depends on 
m, Fd and Rg or Rt and where Adgq —» 0 when Fd 0. In other 
words, W is minimised when the diver is neutrally buoyant at 
the middle-depth of the dive D/2. Out from the minimum, vari¬ 
ations of a few percent of the value of Rt lead to changes of W 
of several orders of magnitude. When the value of R, is rela¬ 
tively far from one, i.e. in cases D and E, W increases linearly 
with it, as expected from Eqs. (O. W is much less sensitive 
to changes of Rg. Variations of several times its value produce 
moderate changes in the value of W. 


An overall view of the locomotion cost dependency on diver 
characteristics is presented in Eigure[3](top) where we plot W as 
a function of both the R, and Rg parameters. Eollowing the grid 
lines corresponding to constant R, values, equivalent to the plot 
in Eigure 12 (bottom), we can observe that the value —» 0 

when /?,—>!. This is connected to the presence of the global 
minimal region of W close to the values Rt - 1 and Rg - 0. 
When Rg ~ 0 and f?, ~ 1, Eq. (fTOl) is satished and case A is re¬ 
placed by case E whose range corresponds exactly to the global 
minimal region of W where W - 2DFd, and no gliding phases 
are present in the dive. A close-up of this global minimal re¬ 
gion is presented in Eigure [3] (bottom) with the cut curve W{Ri) 
corresponding to the constant value Rg = 0.005. A similar be¬ 
haviour, with a minimal region corresponding to case E, can 
also be found for a high drag force (high swim velocity) and/or 
dives to a limited maximal depth, for which Eq. (nnii can also 
be satished. In the extreme condition where Rg = 0, only cases 
D, E and E are possible (no presence of critical depths), and W 
is minimal when 


.. Fd Fd 

1 - <Rt < \ + — 

mg mg 


( 12 ) 


is satished (case E). 

The above results about the minimisation of the thrust work 
W can also be understood qualitatively with a simple reasoning; 
during the dive, thrust has to counterbalance the action of drag, 
which is a non-conservative force, and weight and its coun¬ 
terpart the buoyant force, which are both conservative forces. 
The work spent against the conservative forces in one phase 
(ascent and/or descent) is given back in the opposite phase of 
the dive and is used for gliding. During the glides, the buoy¬ 
ant force increases or decreases progressively with the change 
of depth, produci ng diver acceleration, as observed in penguins 
( Sato et all 2002 ). with a consequent increase in drag. It is pos- 
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Figure 3: Top: Thrust work as function of R[ and Rg for a bottlenose dolphin 
of mass m = 177 Kg and drag constant C = 4.15 kg diving to a maximum 
depth of 100 m with a speed of 1.8 m s“^. Bottom: a constant Rg curve (Rg = 
0.005) close to the minimal W region. 
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sible for the diver to maintain constant velocity, but in this case 
it must increase drag intentionally. In other words, during the 
glides, the potential energy related to gravitation is dissipated 
in the turbulent motion of the water and then lost. To minimise 
this, the glide regions should be reduced as much as possible. 
This means that the critical depths should be kept as close as 
possible to the middle-depth D/2 corresponding to condition of 
neutral buoyancy of the diver around D/2, far from the condi¬ 
tions with a very prolonged glide in one phase of the dive and 
another short glide or no glide at all in the other phase (like in 
case D and E). This is exactly what we found formally from the 
formulae plots presented above. When Rg (or R,) is kept con¬ 
stant, the values (or that minimises W corresponds 
to the critical depths deqD,A - ± and neutral buoyancy 

of the diver at exactly half the maximal depth of the dive D/2. 
Note that this is in fact a well known practice in freediving (a 
human sport). Where both Rg and are varied, the best values 
are ~ 0 and ~ 1, which implies case F with total absence 
of glides and close to neutral buoyancy. 


The previous consideration implies that when Rt value is suf¬ 
ficiently close to one, a decreasing of thrust work can be ob¬ 
tained by exhaling before diving, to reduce the value of Rg. 
This operation also sensitively reduces the amount of available 
stored oxygen, penalising the maximum achievable depth. This 
is not true in the case of phocid seals, where stored oxygen 
is mainly located in the blood and muscular tissue and where 
the lung air repre sents only 5% of available oxygen reserves 
( Koovmani 1985)). Phocid seals are in fact observed to exhale 
before diving ( ^holande3 194(]l: Koovman et al. . 197(]h . This 
practice is generally interpreted as a red uction in risk of decom¬ 


pression sickness v ia alveolar co ll apse (IKoovman et al.L 11972 


FaUce et al.L 119851: iHooker et al.L 120051) . similarly to otariid 


seals, whic h are observed to ex hale in the ascending phase of 


their dives (iHooker et al.L l200-5h . An additional explanation to 


these behaviours could, however, be the intention of reducing 
dive locomotion cost. 


From Figs. |2] and [3] we can see that the measured dolphin’s 
body values = 0.9809 and = 0.0493 considered 
here are close to the minimal region demonstrating dolphins’ 
good adaption for dives around 100 m deep. This good adaption 
is not limited to these depths but also for a larger range value 
as we can see from the plot in Figure 0] Here we plot the value 
of W divided by diver mass and the distance covered (equal to 
2D) as a function of the maximal depth. We see that, except 
for dives close to the surface, the normalised locomotion cost 
is almost constant, demonstrating a good adaption of dolphin 
body characteristics for dives down to 200 m and more. For 
comparison, we plot the normalized locomotion cost for em¬ 
peror and king penguins. In particular we consider the example 
of a specific emperor penguin with m - 30 Kg, = 4 1 and 
C = 0.60 kg m ' (data from Sato et al. ( 2010l) . corresponding 
to the values Rg - 0.1369 and/?, = 1.0069), and a king penguin 
with m = 1 1.5 Kg , Vg - 1.48 1 and C = 0.50 kg m"' (data from 
Sato et al. ( 20021) . corresponding to the values Rg - 0.1322 and 


R, - 1.0069). For all divers we consider the same cruise veloc¬ 
ity of 1.8 m s '. As mentioned above, due to their low average 
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Figure 4: Locomotion cost per unit of covered distance and diver body mass as 
function of the dive deptii for bottlenose dolphins, emperor and king penguins 
swimming at a constant speed of 1.8 m s“'. We consider a bottlenose dolphin 
with mass m = 177 Kg, drag constant C = 4.15 kg m“* and body volumes Rg = 
0.0493 and R, = 0.9809, a king penguin with mass m = 11.4 Kg, C = 0.50 kg 
m“' and body volumes Rg = 0.1322 and R, = 1.0069 and an emperor penguin 
with mass m = 30 Kg, C = 0.60 kg m“* and body volumes Rg = 0.1369 and 
R, = 1.0069. 


tissue density compared to the salty water (/?, > 1), emperor 
penguins are always positively buoyant at any depth and so fall 
into case E, with a glide during the entire ascent. In the case 
of the king penguin, this is not completely true. If the dive is 
too deep, even if they are always positively buoyant, the ani¬ 
mal has to actively swim in the first part of the ascent against 
the drag (c ase C). This is in agreement with the observation 
reported by Sato et al. ( 2002h (see next section for a quantita¬ 
tive comparison between our prediction and these data). As 
expected from Figure|2] for penguins the thrust work is greater, 
and decreases progressively with the depth with a global energy 
economy worse than dolphins. From the consideration above, 
this can be interpreted as consequence of the absence of neu¬ 
tral buoyancy during the dive. Dolphins have neutral buoyancy 
(without considering the drag) at depth c/ 15 m with the pres¬ 
ence of the prolonged glides in both descending and ascending 
phases. In these phases, weight and buoyant forces are used best 
against the drag, resulting in a much lower and the locomotion 
cost. This is not completely true for shallow dives, where the 
prolonged glide during the descent cannot be present and the 
dolphin glides for the entire ascent (case E) causing an increase 
in the dive cost. 

The optimal value minimising W depends only 

marginally on the drag force (on F^l(mg)^ or higher terms). 
This means that dolphin characteristics are not only optimal for 
a wide range of maximal depths but also for a wide range of 
swim velocities. 


3.2. Gliding region extension 

For a defined set of parameters, Eqs. (|7]i can be used to pre¬ 
dict the depths where the diver should begin their glide. We can 
perform a first quantitative comparison for king penguins, for 
which a large set of glide start depth values d'^“‘^' is reported 
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Sato et al.l (120021) . The depth of active swim cessation in the as¬ 


cent depends on pre-dive lung filling and pitch angle assumed 
during the dive. Assuming the maximal air lung filling and 
vertical dives, Eqs. O provide the maximum reachable value 
of (Ia- Considering the kin g penguin charact eristics and their 
typical swim velocity from ISato et al.l (l2002h (m - 11.5 Kg, 
Vg - 1.48 1 and C = 0.50 kg m“', v = 2 m s“'), we find that cIa 
values should be less than 111m. This value is well in agree- 
ment with the ob served which are all less than 100 m 

( Sato et al. . 2002). 


A much more quantitative comparison can be done with dol¬ 
phins for which observations of pure vertical dives are avail¬ 
able. For a 100 m deep dive by the same dolphin considered 
in the previous section, Skrovan et. al observed a glide during 
the descent starting at a depth of = 67.5 m preceded by 
an active swimming phase with a cruise velocity of 1.7 m s“', 
and a glide during the ascent starting at = 5.5 m preceded 
by an active swimming phase with a cruise velocity of 1.9 m 
s *. For the same conditions, our model predicts the values 
do - 30.5 m and = 8 m, which are significantly different 
from the measured values. A possible partial pre-dive filling of 
the lungs, and then a reduced value of Rg, causes a reduction in 
both critical depths, ruling this out as an explanation for the in¬ 
consistency with the observed values. The most probable cause 
of this discrepancy is an incorrect evaluation of the drag force, 
i.e., the value of the drag constant C or of the ratio A between 
the active and passive swimming drag force. This can be eas¬ 
ily understood from Fqs. (0. An increasing of causes an 
increasing of do and a decreasing of dA- Other parameters of 
Fq. 0 cannot be varied easily. 


For our predictions, we use the drag constant for uninstru¬ 
mented dolphins (C = 4.15 kg m '), but the observations of the 
gliding region s are re lative to instrumented dolphins for which 


Skrovan et al.l (Il999h observe a drag increase of more than 4 


times (corres ponding to C = 16.9 kg m *, extracted from the 
data points in Skrovan et al. ( 1999l) f. With this drag constant, 
we should observe almost no gliding phases except for the very 
end of the ascent. The same results can be obtained by tak¬ 
ing a value of A several times larger than one. The problem 
is, in fact, more complex because carrying instruments could 
aff ect the animal behaviour i n an non-trivial way, as discussed 
by van der Hoop et all ( 2014 ). In addition in our simple model 


we do not consider the stroke-and-glide swimming that could 
introduce an additional dependency on v of the average values 
of C and A. A deeper understanding of passive and active drag 
is required for a better prediction of glide phase extent. 


We can, however, provide the maximal value that C and A 
can assume from comparison between obverved values and the 
simple formulae of the critical depths (Fqs. 0). Their single 
values cannot however be disentangled from each other. In¬ 
dependently from the possible approximations, we find that A 
cannot be higher than 1.5 (i.e. C cannot be 1.5 time larger 
than the value for uninstrumented dolphins), much less than the 
value considered in some past studies which indicates 4 3-7 


(Lighthilll 1971; 

Webbl 

1975 

Fish, 

1998; Skroyan et al.. 1999; 

Anderson et al., 2001 

Fish et all 2014. 


3.3. Swim velocity and optimisation 


As announced in Sec. 12.41 the choice of the swim velocity 
is a crucial point for the total energy cost E due to the basal 
metabolism, which contribution is proportional to the total div¬ 
ing time T = 2D/v, and the drag, whose increases with the 
square of the swim velocity v. 

In the general case we have a non-trivial relationship between 
E, given by Fq. (fTTb . and v. This is not completely true in cases 
D, F and F, where a simpler expression for E(D, v) can be found 
due to absence of any critical depth in the expression of W. For 
cases D and F the diver is always negatively or positively buoy¬ 
ant and the expression of W is simple (Fqs. 0). The minimi¬ 
sation of E with respect to v leads to the optimal cruise velocity 



which is independent of maximal d i ve dep th. This expression 
was previously found bv ISato et al.l (120101) for the study of the 
optimal speed of emperor penguins, which are always posi¬ 
tively buoyant and are included in the case F where Fi, - F„ > 

Ff- 

For case F, the expression of the optimal speed is slightly 
different due to the thrust work required for both descent and 
ascent. Considering Fq. (fTTT l and the expression of W for case 
F, Eq. 0, we obtain an optimal velocity solution very similar 
to Eq. (foi l with the denominator term C in Eq. (fT3l l replaced 
by 2C. N.B. this is also the optimal velocity for a horizontal 
displacement. 

No simple analytical formula of can be obtained for the 
other cas es and numerical methods (Mathematica software in 
our case d Wolfram Research . 2015 )1 must be applied for solv¬ 
ing the equation dEjdv - 0. To better compare our results with 
previous studies, instead of E we consider the cost of trans¬ 
port per units of mass and covered d istance COT - EHID m). 
In addition to m,R,,Rg and C from Skrovan et al. ( 19991) . the 
metabolic and propuls ion efficiencies = 0.25 and Sp = 0.86 
(IFish and Rohri Il999h are used and where aerobic efficiency 
only is taken into account for s. This is ju stified by post¬ 
dive b lood lactate measurements in dolphins ( Williams et al.l 
19991) showing an increase in anaerobic process only after sig¬ 


nificantly long dives (longer than 200 s), not considered here. 

For the value of B, we co nsider the basal metabolic rate re¬ 
ported by Yazdi et ^ ( 1999 ) {B/m - 2.15 W Kg“*) measured 
for trained dolphins in a pool but considering some precautions. 
During the dives, dolphins, like other breath-hold divers, are 
subject in fact to a physiological response to diving (diving re¬ 
flex) with bradycardia, reduction and redistribution of the blood 
flow, etc. Specific ally, a reduction in hea rt rate to 63.4% is ob¬ 
served in dolphins dWilliams et al.ll 19991) with an expected sub¬ 
stantial reduction of B. To roughly take into account the diving 
reflex, we also examine the expected value of COT correspond¬ 
ing to half of the basal metabolic rate value reported in litera¬ 
ture (B/m = 1.07 W/Kg), which is also approxi mately equal to 


basal metabolism predicted by the Kleiber’s law (lKleiberlll947 
Yazdi et al.Lll999l) . 
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Figure 5: Top: Cost of transport as a function of dive velocity v for a bot- 
tlenose dolphin of mass m = 177 Kg, drag constant C = 4.15 kg and 
body volumes Rg = 0.0493 and Rr = 0.9809 diving to a maximum depth of 
D = 100 m. Different curves con'espond to different assumptions for the basal 
metabolic rate B and the active to passive drag ratio A. The presence of critical 
depths depends on the value of v (both critical depths present: case A; only 
one present in the ascent: case C; none present: case F). Bottom: Close-up of 
the minimum region where the dependency on D is more visible and where we 
assumed Bjm = 1.07 W kg“^ and A= \. 


In Figure |5] we plot the COT for the same dolphin individ¬ 
ual considered in the previous section for a 100 m depth dive 
for different combinations of values of Bjm (Bjin - 2.15 and 
1.07 W Kg“’) and A (1 and 1.5). Changes in dive velocity de¬ 
termine the existence or not of critical depths and so which dive 
case applies. Here in particular, cases A, C and F must be con¬ 
sidered. For Bjm - 2.15 W Kg ' and d = 1 the optimal ve¬ 
locity, solution of the equation dEjdv - 0, is = 2.3 m 
s '. This valu e is quite different fro m the value of about 1.8 m 


s 'reported by Skrovan et al. ( 19991) . When the basal metabolic 
rate value of Bjm - 1.07 W Kg ' is considered to include the 
diving response, an optimal velocity of = 1.95 m s“' is 
obtained, very close to the measured values. An underestima¬ 
tion of the drag force can also produce an artificial decreasing 
of We remember that we consider here a value of the drag 
constant for uninstrumented dolphins when the observations are 
made on instrumented dolp hins. If we take the v alue of C for 
instrumented dolphins from (ISkrovan et al.lll999h fC = 16.9 kg 
m '), we get the very low value of v*®*' = 1.3 m s“' - very dif¬ 
ferent from the observed value, as we found for critical depths 
(Section I3.2l i. Similarly, if we consider larger values of A, a 
smaller is obtained. In fact, when the A upper limit of 1.5 
is adopted (obtained from the observation on the critical depths 
in Sec. 13.21 1. the value = 1.95 m s 'is found (see Fig. |5]l, 
the same value obtained using half the B value from the litera¬ 
ture. 

The cause of the disagreement between observed and pre¬ 
dicted swim velocity when the literature value of B is used is 
probably due to a combination of the effect of the diving re¬ 
sponse and the underestimation of drag. However, considering 
i) the recent studies that indicate a s wim speed reduction of onl y 


11% for instrumented dolphins (Ivan der Hoop et al 


(carrying diff e rent a pparatus than in IWilliams et al. 


2014), 


19991) : 


in some recent articles ( 

Hind and Gurnev, 1997; Barrett et al.l 

19991 

Boraziani and Sotironoulosl 2008 

) and iii) the observed 

heart rate reduction of diving dolphins ( 

Williams et al.. 1999t). 


we suggest that the diving reflex plays the major role and a re¬ 
duced value of B should be considered in favour of an increase 
of /I or C values. For simplicity, the following calculations use 
the values B/m = 1.07 W kg“', A - I and C = 4.15 kg m“' 
(from uninstrumented dolphins). 

The dependency of the COT on maximal dive depth is 
marginal, as it can be observed in Figure |5] (lower panel) where 
we plot the COT for dives with D - 25,50,100 m. Optimal 
dive velocity dependency on maximal dive depth and the cor¬ 
responding value of COT are presented in Figure Different 
depths correspond to different cases. As expected, large val¬ 
ues of D fall into case A. For the particular dolphin individual 
considered, dives of less than 33 m fall into case C (no pro¬ 
longed glide in the descent). Presence or absence of prolonged 
glides and dependency on their extent produce a non-trivial de¬ 
pendency of the optimal speed on the maximal depth. 
varies about 10%, in the 1.8 and 2.0 m s ' range, similar to 
the observed values ( Skrovan et al. . 1999h . Differently to W 


(Fig. |4]i, the cost of transportation calculated for v = has a 
minimum with respect to the maximal dive depth D, which in 
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Figure 6: Optimal velocity dependency (top) and coiTespondent COT (bottom) 
on the maximal dive depth for a bottlenose dolphin of mass m = 177 Kg, drag 
constant C = 4.15 kg m“* and body volumes Rg = 0.0493 and R, = 0.9809. 
Solution case A or C has to be considered. Both cases are represented and the 
pertinent values of which depends on D, are indicated by thick lines. 


this case is around D ^ AQm. 

We can study the dependency of and the corresponding 
COT on changes of Rt and Rg ratios, as we did for locomotion 
work. In this case, however, we have a more complex problem. 
^hest j-jjg solution of dEjdv - 0, where E must correspond to 
the pertinent case (A to F), whose selection depends on the Ed 
value and then on v. 






Figure 7: Optimal velocity dependency (top) and correspondent COT (bottom) 
on the ratio R[ for a bottlenose dolphin of mass m = 177 Kg and drag constant 
C = 4.15 kg m“^ diving to a maximum depth of D = 100 m keeping Rg constant 
at 0.0493 (= 


As we can see in Figure|2](top panel), when we vary Rt keep¬ 
ing Rg = constant, varies up to 15% in the region 

Rt ~ I (when prolonged glides are present). When Rt is rela¬ 
tively far from one, the diver is extremely positively or nega¬ 
tively buoyant (cases D and E) and the optimal velocity expres¬ 
sion is given by the Eq. (fTST l. independent of the magnitude of 
Rt and Rg values. 

As we can see in Eigure|3(top panel), variation in Rg keep¬ 
ing Rt constant causes the switch from one case to another with 
extreme cases E and D or E, depending on the considered value 
of Rt- Close to neutral buoyancy {Rt values enough close to 
one), case E in fact replaces case D. This is visible also in Eig- 
ure|9](top panel) where the variation of with respect on Rt 
is presented for the extreme case Rg = 0. 

The optimal cost of transportation corresponding to v = 
for different values of Rt and Rg is presented in the lower pan- 
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Figure 8: Optimal velocity dependency (top) and coiTespondent COT (bottom) 
on the ratio Rg for a bottlenose dolphin of mass m = 111 Kg and drag constant 
C = 4.15 kg m“^ diving to a maximum depth of D = 100 m keeping R{ constant 
at 0.9809 (= ) and 0.9950. 


Figure 9: Optimal velocity dependency (top) and con'espondent COT (bottom) 
on the ratio R, for a bottlenose dolphin diving to a maximum depth of D = 
100 m for the extreme case Rg = 0. 
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Figure 10: Breath-hold diver velocities for a dive 100 m depth as function of 
diver mass. Theoretical predictions with different assumpti ons are compared to 
available observations in mammals iWatanabe et al.Ll201lh . 


els of Figures|7]-|9] Dependency of the COT on ratios /?, and 
Rg is similar to W, with a global minimum close to the region 
Rg - Q,Rt = I, i.e. for almost neutrally buoyant divers with 
a very small gas-filled body volume. This confirms the advan¬ 
tage of exhaling before diving dive, observe d in phocid seals . 


whose Rt - pjp t value can be close to one (lAoki et al.L 1201 1 


Sato et al. . 2013h . We remember that unlike other marine mam¬ 


mals, air stored in phocid seal lungs represents only the 5% of 
the total stored oxygen, which is mainly located in the blood 
and in muscular tissue. In this case, air exhalation does not 
reduce too much the total available oxygen but does allow for 
reducing the COT. 


As for W, the minimal region corresponds to the case F, for 
which the expression of COT becomes very simple and does 
not depend on physical characteristics of the diver; 


COT^iv = 

m V 4s 


(14) 


For the extreme condition where Rg - 0, unlike W, in ad¬ 
dition to the possible cases D, E and F, we have two transition 
regions corresponding to case F limits with cases D or E, where 
= 0 or (fo = D, respectively. Unlike W, we are in fact adjust¬ 
ing V to minimise the COT for each value of Rt which influences 
case selection. 


The presence of the COT mi nimum close to Rt - 1 when 
R„ - 0 confirms the results from Miller et al. ( 2Q12h : Sato et al 


(120 13h . where, when Rg - 0 (equivalent to considering buoy¬ 
ancy not dependent on d), the resulting locomotion cost is min¬ 
imal at neutral buoyancy but where no variation of the buoyant 
force is considered. Here we extend this result to a larger range 
of values of Rt corresponding to the region defined by Eq. (fTSli 
with 


Fd = = 


Ce^B^ 


(15) 


3.4. Mass dependency of optimal dive velocity 

A more general comparison between the prediction of our 
model and observed va lues can be don e by co nsidering differ¬ 
ent diver masses. Like I Watanabe et alJ (1201 11) . we introduce a 
mass dependency on the different parameters of Eqs. (|4]-[TTJ. 
Eor the drag constant C we consider a variation of from 
the dependency of the diver surface and the drag co¬ 

efficient (m“^^'^). Eor the basal metabolic rate we take the 
classic dependency of from the Kleiber’s law, keeping 
as reference half the value of B (B relative = 1/2) of rest¬ 
ing dolphins at the surface ( Yazdi et al. . 1999ll to take into ac- 
count the physiological respons e to diving. Erom observations 
( Koovman and Ponganis . 1998h and basic principles, it is dif¬ 
ficult to make an estimation of lung volume variation with the 
mass. We consider then, no dependency on m for Rg as for the 
efficiency s and we retain the value 4=1. 

In Eigure [TOl we plot for a dive with D — 100 m as 
function of m for a large range of diver masses, all correspond¬ 
ing to case A w here both critica l depth s are present. As we can 
see, contrary to lWatanabe et all (1201 ill , no exponential law can 
be extrapolated. Com paring our predictions t o the observed ve¬ 
locities of mammals ( Watanabe et al. . 201 ih (which generally 
glide in both ascending and descending phases) we see that our 
predictions are in quite good agreement for small values of m 
but are systema t ically higher for larger masses. In their study 


Watanabe et al. ( 201ll) suggest a possible reason for the ob¬ 


served large difference in seabird and mammal swim velocities 
(both endotherms) could be due to the large difference in ther¬ 
moregulatory costs between the species. Similarly, a reduction 
in thermic loss coupled to a strong diving reflex could explain 
large mammals’ small swim velocities. The discrepancy be¬ 
tween predicted and observed velocities could also come from 
the approximations made in our model, such as the assumptions 
made for A and B. When a greater value of 4 (4 = 1.5 instead 
of 1) and any reduction of B due to the diving reflex {B relative 
= 1) are considered, the disagreement between predictions and 
observations increases. This suggests that the origin of this dis¬ 
crepancy is elsewhere, such as i) the value of B discussed above, 
or ii) neglect of possible accelerations during prolonged glides 
or iii) rough approximations of the dependency of the differ¬ 
ent parameters, in particular C and Rg, on diver mass. Eurther 
studies are required to clarify this point. 


4. Conclusions 


We presented a theoretical model for breath-hold animal 


dives. In accordance to previous studies (Wilson et al.l 1992 

Skrov 

anetall 1999: Hansen and Ricklef 

2004 

; ISato et al. 

201C; 

AokietakL 2011; Watanabe et al.. 

ZOlll 

Miller et al. 

2012) 

, starting from basic principles of forces acting on the 


diver, we calculated the work of the thrust to evaluate ener¬ 
getic locomotion dive cost. Unlikely to previous studies, , we 
were able to take into account accurately for the first time the 
presence of prolonged glides during the ascent and descent, cal¬ 
culate their extent and study their impact on the total locomo¬ 
tion cost. In particular we analyse locomotion cost dependency 


12 

























































































on different dive parameters, namely diver body characteris¬ 
tics (mass, lung volume, buoyancy, etc.) and dive character¬ 
istics (cruise velocity and maximum depth). To calculate total 
dive energy cost, we also include the basal metabolism and ef¬ 
ficiency for transforming chemical energy to propulsion. We 
studied total dive cost dependency on cruise velocity and the 
choice of optimal swim speed to minimise dive cost as function 
of the different parameters. 


We demonstrate that both locomotion and total energy cost 
are minimised for divers which pass through a condition of neu¬ 
tral buoyancy during the dive, generally implying the presence 
of prolonged gliding phases in both ascent and descent. This 
is in a greement with past results ( Miller et al. . 2012t Sato et al 


2013h where, when the buoyant force is considered constant 


during the dive, energy cost is minimised for neutral buoyancy. 
In particular, our model confirms good physical adaption of dol¬ 
phins for dives of 10 - 200 m deep due to their small body 
tissue density over mass ratio, contrary to the case of emperor 
king penguins. The presence of prolonged glides implies a non¬ 
trivial dependency of optim al speed on maximal dive depth, and 
exten ds previous findings ( Sato et all 2010l : Watanabe et al 


201 ll) who found no dependency on dive depth. Locomotion 


and total energy cost of the dive are further diminished by re¬ 
ducing gas-filled body parts volume for divers close to neutral 
buoyancy. In particular, this provides an additional possible ex¬ 
planation of the pre-diving exhalation of air observed in pho- 
cid seals to minimise dive energy cost, rather than the current 
explanation from the literature of decompression sickness risk 
reduction. 


We also successfully compare the extent of gliding phases 
with observations of penguin and dolphin dives. In particular, 
from the quantitative comparison b etween the da t a for v erti- 
cal dives of a bottlenose dolphin dSkrovan et al.L Il999h and 
our model predictions, we estimate an upper limit of the ra¬ 
tio of active to passive drag of 1.5, confir ming the find i ngs of 
studies with robotic fish and sim ulations ( Barrett et al. . 19991 : 
Boraziani and Sotiropoulosl 2008 ). 


The comparison between our findings and observed swim ve¬ 
locities of breath-hold mammals with different masses suggests 
that the physiological response to diving could play an impor¬ 
tant role in the choice of cruise velocity but it also shows the 
limitations of our model and its approximations. In particular, 
the assumption of a constant cruise velocity during the dive and 
the dependency of the drag on the swim speed, including a bet¬ 
ter involvement of the stroke-and-glide swimming mode, and 
on other parameters, should be revisited and investigated. 
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